2. Free Fields 



"The career of a young theoretical physicist consists of treating the harmonic 

oscillator in ever-increasing levels of abstraction." a . , „ , 

bianey Coleman 

2.1 Canonical Quantization 

In quantum mechanics, canonical quantization is a recipe that takes us from the Hamil- 
tonian formalism of classical dynamics to the quantum theory. The recipe tells us to 
take the generalized coordinates q a and their conjugate momenta p a and promote them 
to operators. The Poisson bracket structure of classical mechanics morphs into the 
structure of commutation relations between operators, so that, in units with h — 1, 

[qa,q b ] = [p a ,p b ]=0 

[<la,p b }=i6 b a (2.1) 

In field theory we do the same, now for the field (f) a {x) and its momentum conjugate 
it b (x). Thus a quantum field is an operator valued function of space obeying the com- 
mutation relations 

iMx)A b m = i* a (x),Ay)] = o 

[Mx),Ay)]=^ 3 \x-y)5 b a (2.2) 

Note that we've lost all track of Lorentz invariance since we have separated space x 
and time t. We are working in the Schrodinger picture so that the operators <f) a (x) and 
ir a (x) do not depend on time at all — only on space. All time dependence sits in the 
states which evolve by the usual Schrodinger equation 

(2.3) 

We aren't doing anything different from usual quantum mechanics; we're merely apply- 
ing the old formalism to fields. Be warned however that the notation for the state 
is deceptively simple: if you were to write the wavefunction in quantum field theory, it 
would be a functional, that is a function of every possible configuration of the field 0. 

The typical information we want to know about a quantum theory is the spectrum of 
the Hamiltonian H . In quantum field theories, this is usually very hard. One reason for 
this is that we have an infinite number of degrees of freedom — at least one for every 
point x in space. However, for certain theories — known as free theories — we can find 
a way to write the dynamics such that each degree of freedom evolves independently 
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from all the others. Free field theories typically have Lagrangians which are quadratic 
in the fields, so that the equations of motion are linear. For example, the simplest 
relativistic free theory is the classical Klein-Gordon (KG) equation for a real scalar 
field (f>(x, t), 

d^(f) + m 2 (f) = (2.4) 

To exhibit the coordinates in which the degrees of freedom decouple from each other, 
we need only take the Fourier transform, 



Then (f>(p, t) satisfies 



q2 

_ + (^ + m2) ) 0(£ t ) = O (2.6) 



Thus, for each value of p, 0(p, t) solves the equation of a harmonic oscillator vibrating 
at frequency 



LO 



ff =+Vp 2 + ™ 2 (2-7) 



We learn that the most general solution to the KG equation is a linear superposition of 
simple harmonic oscillators, each vibrating at a different frequency with a different am- 
plitude. To quantize <j)(x, t) we must simply quantize this infinite number of harmonic 
oscillators. Let's recall how to do this. 

2.1.1 The Simple Harmonic Oscillator 

Consider the quantum mechanical Hamiltonian 

H = \p 2 + \uj 2 q 2 (2.8) 

with the canonical commutation relations [q,p] = %■ To find the spectrum we define 
the creation and annihilation operators (also known as raising/ lowering operators, or 
sometimes ladder operators) 

lo % | [uj i 

2 q + 7S p ' a V2 q -7S P (2 ' 9) 



which can be easily inverted to give 



! a + a f ) , p=-i A /^(a-a t ) (2.10) 



2a; v ' 1 1 V 2 
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Substituting into the above expressions we find 

[a, a f ] = 1 

while the Hamiltonian is given by 

H = ^u(aa) + at a) 



(2.11) 



(2.12) 



One can easily confirm that the commutators between the Hamiltonian and the creation 
and annihilation operators are given by 



[H, at 



ooa 



f and [H, , 



-ooa 



(2.13) 



These relations ensure that a and take us between energy eigenstates. Let \E) be 
an eigenstate with energy E, so that H\E) = E\E). Then we can construct more 
eigenstates by acting with a and at , 



Hat \E) = (E + oo)at \E) , Ha \E) = (E — oo)a \E) 
So we find that the system has a ladder of states with energies 

. . . , E - u, E, E + u, E + 2uj, . . . 



(2.14) 



(2.15) 



If the energy is bounded below, there must be a ground state |0) which satisfies a |0) = 0. 
This has ground state energy (also known as zero point energy), 

ff|0> = |a;|0> (2.16) 

Excited states then arise from repeated application of a\ 

|n) = (a t ) n |0) with H \n) = (n + \)u \n) (2.17) 

where I've ignored the normalization of these states so, (n\ n) ^ 1. 

2.2 The Free Scalar Field 

We now apply the quantization of the harmonic oscillator to the free scalar field. We 
write and n as a linear sum of an infinite number of creation and annihilation oper- 
ators at and a$, indexed by the 3-momentum p, 



0(f) 
tt(x) 



d 3 p 



(2tt)3 ^ 
d 3 p 

(2^p ( ^ 



a^e 



'" x + ale~ if * 




a*e if - s - aie- if - s 



(2.18) 
(2.19) 
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Claim: The commutation relations for <p and ir are equivalent to the following com- 
mutation relations for ap and at 

[0(f), <P{y)\ = Hf), <y)\ = [a-, a,-] = [aj, aj] = 



Proof: We'll show this just one way. Assume that [a-, at] = (2ir) 3 5^(p — q). Then 



(2tt) 6 2 Y up 

I" d 3 p (-t) , ip-(x-y) _ Jp-(y-x) 

" J (2vr)3 2 [ 6 



) 



(2.21) 
□ 



The Hamiltonian 

Let's now compute the Hamiltonian in terms of a- and at. We have 



H=\j#x 7r 2 + (V0) 2 + m 2 2 



d 3 x d 3 p d 3 g 



1 f d i x d 6 i 
= 2 J (2^ 



{ipa-jv- s - ipate-**) ■ (iqa-e 1 ^ - if ale* 1 *) 



)(a-e^-ate-^) 



I [ d 3 p 1 



(2vr)3 Wp 



a f e** + ate-^)(a ? -.e^ + ate"^) 

+ p 2 + m 2 ) (a^a^ + atal-) + (cj 2 . + p 2 + m 2 ) (a^at + ata p ->) 



where in the second line we've used the expressions for and 7r given in (2.18) and 
(2.19); to get to the third line we've integrated over d 3 x to get delta-functions 5^(p±q) 
which, in turn, allow us to perform the d 3 q integral. Now using the expression for the 
frequency uj'i = p 2 + m 2 , the first term vanishes and we're left with 



H = 



d 3 p 

(2^ 
d 3 p 



a p a l +a l a p 



(27T) 



3 U ? 



ala,+ ^n) 3 6^ 3 \0) 



(2.22) 
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Hmmmm. We've found a delta-function, evaluated at zero where it has its infinite 
spike. Moreover, the integral over uif diverges at large p. What to do? Let's start by 
looking at the ground state where this infinity first becomes apparent. 

2.3 The Vacuum 

Following our procedure for the harmonic oscillator, let's define the vacuum |0) by 
insisting that it is annihilated by all a-, 

a p \0) = Vp (2.23) 

With this definition, the energy E of the ground state comes from the second term in 
(2.22), 



H\0) =E \0) 



d 3 p\u^\0) 



0) = oo |0) (2.24) 



The subject of quantum field theory is rife with infinities. Each tells us something 
important, usually that we're doing something wrong, or asking the wrong question. 
Let's take some time to explore where this infinity comes from and how we should deal 
with it. 

In fact there are two different oo's lurking in the expression (2.24). The first arises 
because space is infinitely large. (Infinities of this type are often referred to as infra-red 
divergences although in this case the oo is so simple that it barely deserves this name). 
To extract out this infinity, let's consider putting the theory in a box with sides of 
length L. We impose periodic boundary conditions on the field. Then, taking the limit 
where L — > oo, we get 

fL/2 j-L/2 

(2tt) 3 5 (3) (0) = lim / d 3 x e^L = lim / d 3 x = V (2.25) 

L ^°°J-L/2 V L ^°°J-L/2 

where V is the volume of the box. So the 5(0) divergence arises because we're computing 
the total energy, rather than the energy density So- To find So we can simply divide by 
the volume, 



Eo f d 3 p 1 



(2.26) 



which is still infinite. We recognize it as the sum of ground state energies for each 
harmonic oscillator. But S — > oo due to the \p\ — > oo limit of the integral. This is 
a high frequency — or short distance — infinity known as an ultra-violet divergence. 
This divergence arises because of our hubris. We've assumed that our theory is valid 
to arbitrarily short distance scales, corresponding to arbitrarily high energies. This is 
clearly absurd. The integral should be cut-off at high momentum in order to reflect 
the fact that our theory is likely to break down in some way. 
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We can deal with the infinity in (2.24) in a more practical way. In physics we're 
only interested in energy differences. There's no way to measure Eq directly so we can 
simply redefine the Hamiltonian by subtracting off this infinity, 

so that, with this new definition, H |0) = 0. In fact, the difference between this Hamilto- 
nian and the previous one is merely an ordering ambiguity in moving from the classical 
theory to the quantum theory. For example, if we defined the Hamiltonian of the har- 
monic oscillator to be H = (1 /2)(ujq — ip){ujq + ip) , which is classically the same as our 
original choice, then upon quantization it would naturally give H = ua^a as in (2.27). 
This type of ordering ambiguity arises a lot in field theories. We'll come across a number 
of ways of dealing with it. The method that we've used above is called normal ordering. 

Definition: We write the normal ordered string of operators . . . (p n (x n ) as 

:0i(£i)...0„(£„) : (2.28) 

It is defined to be the usual product with all annihilation operators a- placed to the 
right. So, for the Hamiltonian, we could write (2.27) as 

■H:= J -^Ujta^ap (2.29) 

In the remainder of this section, we will normal order all operators in this manner. 
2.3.1 The Cosmological Constant 

Above I wrote "there's no way to measure E directly" . There is a BIG caveat here: 
gravity is supposed to see everything! The sum of all the zero point energies should 
contribute to the stress-energy tensor that appears on the right-hand side of Einstein's 
equations. We expect them to appear as a cosmological constant A = E /V, 

R^u - = -BnGT^ + A 9lu , (2.30) 

Current observation suggests that 70% of the energy density in the universe has the 
properties of a cosmological constant with A ~ (10~ 3 eV) 4 . This is much smaller than 
other scales in particle physics. In particular, the Standard Model is valid at least up 
to 10 12 eV. Why don't the zero point energies of these fields contribute to A? Or, if 
they do, what cancels them to such high accuracy? This is the cosmological constant 
problem. No one knows the answer! 
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2.3.2 The Casimir Effect 



"I mentioned my results to Niels Bohr, during a walk. That is nice, he 
said, that is something new... and he mumbled something about zero-point 



energy. 



Hendrik Casimir 



Using the normal ordering prescription we can happily set E = 0, while chanting 
the mantra that only energy differences can be measured. But we should be careful, for 
there is a situation where differences in the energy of vacuum fluctuations themselves 
can be measured. 

To regulate the infra-red divergences, we'll make the x l direction periodic, with size 
L, and impose periodic boundary conditions such that 



(x + Ln) 



(2.31) 



with n = (1,0,0). We'll leave y and z alone, but remember 
that we should compute all physical quantities per unit area 
A. We insert two reflecting plates, separated by a distance 
(I < L in the x 1 direction. The plates are such that they 
impose 0(x) = at the position of the plates. The presence of 
these plates affects the Fourier decomposition of the field and, 
in particular, means that the momentum of the field inside the 
plates is quantized as 



V 



~d 



neZ A 



(2.32) 



Figure 3: 



For a massless scalar field, the ground state energy between the plates is 
E(d) 



E 



dpydp z 1 / fWK \ 



+ Py+Pt 



A (2tt) 2 2 V V d J 

while the energy outside the plates is E(L — d). The total energy is therefore 

E = E(d) + E(L - d) 



(2.33) 



(2.34) 



which - at least naively - depends on d. If this naive guess is true, it would mean 
that there is a force on the plates due to the fluctuations of the vacuum. This is the 
Casimir force, first predicted in 1948 and observed 10 years later. In the real world, the 
effect is due to the vacuum fluctuations of the electromagnetic field, with the boundary 
conditions imposed by conducting plates. Here we model this effect with a scalar. 
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But there's a problem. E is infinite! What to do? The problem comes from the 
arbitrarily high momentum modes. We could regulate this in a number of different 
ways. Physically one could argue that any real plate cannot reflect waves of arbitrarily 
high frequency: at some point, things begin to leak. Mathematically we want to find 
a way to neglect modes of momentum p ^> a^ 1 for some distance scale a <^ d, known 
as the ultra-violet (UV) cut-off. One way to do this is to change the integral (2.33) to, 




ri7T\ 2 



+P 2 y +P 2 z ) e~ a ^M 2+p y +p * (2.35) 

which has the property that as a — > 0, we regain the full, infinite, expression (2.33). 
However (2.35) is finite, and gives us something we can easily work with. Of course, 
we made it finite in a rather ad-hoc manner and we better make sure that any physical 
quantity we calculate doesn't depend on the UV cut-off a, otherwise it's not something 
we can really trust. 

The integral (2.35) is do-able, but a little complicated. It's a lot simpler if we look 
at the problem in d — 1 + 1 dimensions, rather than d — 3 + 1 dimensions. We'll find 
that all the same physics is at play. Now the energy is given by 



oo 



n=l 

We now regulate this sum by introducing the UV cutoff a introduced above. This 
renders the expression finite, allowing us to start manipulating it thus, 



oo 



n=l 



1 d \ ^ -an-K Id 

2da ^ 

71=1 

1 d I 



2 da 1 - e- a */ d 
7i e an/d 



2d (e a */ d - I) 2 

= + < 2 - 37 > 

where, in the last line, we've used the fact that a <^ d. We can now compute the full 
energy, 

E M = E 1+lW + E 1+1 (L - d) = ^ - i (i + jLj) + 0(a>) (2.38) 
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This is still infinite in the limit a — > 0, which is to be expected. However, the force is 
given by 

dEi+i n 



dd 24d 2 



+ ... (2.39) 



where the . . . include terms of size d/L and a/d. The key point is that as we remove 
both the regulators, and take a — > and L — > oo, the force between the plates remains 
finite. This is the Casimir force 2 . 

If we ploughed through the analogous calculation in d — 3 + 1 dimensions, and 
performed the integral (2.35), we would find the result 

\3E_ 7T 2 

A~dd ~ 480# ( 0) 

The true Casimir force is twice as large as this, due to the two polarization states of 
the photon. 

2.4 Particles 

Having dealt with the vacuum, we can now turn to the excitations of the field. It's 
easy to verify that 

[if, at] = upal and [H, a f ] = -u p a^ (2.41) 

which means that, just as for the harmonic oscillator, we can construct energy eigen- 
states by acting on the vacuum |0) with at. Let 

i£> = 4i°> ( 2 - 42 ) 

This state has energy 

H\p)=wp \p) with up = p 2 + m 2 (2.43) 

But we recognize this as the relativistic dispersion relation for a particle of mass m and 
3-momentum p, 

El = p 2 + m 2 (2.44) 



2 The number 24 that appears in the denominator of the one-dimensional Casimir force plays a more 
famous role in string theory: the same calculation in that context is the reason the bosonic string lives 
in 26 = 24 + 2 spacetime dimensions. (The +2 comes from the fact the string itself is extended in one 
space and one time dimension). You will need to attend next term's "String Theory" course to see 
what on earth this has to do with the Casimir force. 
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We interpret the state \p) as the momentum eigenstate of a single particle of mass m. 
To stress this, from now on we'll write Ep everywhere instead of up. Let's check this 
particle interpretation by studying the other quantum numbers of \p). We may take 
the classical total momentum P given in (1.46) and turn it into an operator. After 
normal ordering, it becomes 

P = - J d 3 xnV(f> = J ^ v 4S~ (2-45) 

Acting on our state \p) with P, we learn that it is indeed an eigenstate, 

P\p)=p\p) (2.46) 

telling us that the state \p) has momentum p. Another property of \p) that we can 
study is its angular momentum. Once again, we may take the classical expression for 
the total angular momentum of the field (1.55) and turn it into an operator, 

r = e ilk [ d 3 x (J°y k (2.47) 



It's not hard to show that acting on the one-particle state with zero momentum, 
J 1 \p = 0) = 0, which we interpret as telling us that the particle carries no internal 
angular momentum. In other words, quantizing a scalar field gives rise to a spin 
particle. 

Multi-Particle States, Bosonic Statistics and Fock Space 

We can create multi-particle states by acting multiple times with a''s. We interpret 
the state in which n a^'s act on the vacuum as an n-particle state, 



\Pl,---,Pn) 



4-... a*, |0> (2.48) 



Because all the gt's commute among themselves, the state is symmetric under exchange 
of any two particles. For example, 

\p,q) = \q,p) (2.49) 

This means that the particles are bosons. 

The full Hilbert space of our theory is spanned by acting on the vacuum with all 
possible combinations of cr's, 

|o> , 4i°) » 44 1°) > 444 io>... (2.50) 
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This space is known as a Fock space. The Fock space is simply the sum of the n-particle 
Hilbert spaces, for all n > 0. There is a useful operator which counts the number of 
particles in a given state in the Fock space. It is called the number operator N 



and satisfies N \pi, . . . ,p n ) = n \pi, . . . ,p n ). The number operator commutes with the 
Hamiltonian, [N, H] = 0, ensuring that particle number is conserved. This means that 
we can place ourselves in the n-particle sector, and stay there. This is a property of 
free theories, but will no longer be true when we consider interactions: interactions 
create and destroy particles, taking us between the different sectors in the Fock space. 

Operator Valued Distributions 

Although we're referring to the states \p) as "particles", they're not localized in space 
in any way — they are momentum eigenstates. Recall that in quantum mechanics the 
position and momentum eigenstates are not good elements of the Hilbert space since 
they are not normalizable (they normalize to delta-functions). Similarly, in quantum 
field theory neither the operators (p(x), nor are good operators acting on the Fock 
space. This is because they don't produce normalizable states. For example, 



(0| a p al |0) = (p\p) = (2tt) 3 5(0) and (0\<j>(x) <p{x) |0) = (x \ x) = 5(0) (2.52) 



They are operator valued distributions, rather than functions. This means that al- 
though <f)(x) has a well defined vacuum expectation value, (O|0(af) |0) = 0, the fluc- 
tuations of the operator at a fixed point are infinite, (0| <f)(x)<p(x) |0) = oo. We can 
construct well defined operators by smearing these distributions over space. For exam- 
ple, we can create a wavepacket 



which is partially localized in both position and momentum space. (A typical state 
might be described by the Gaussian (p(p) = exp(— p 2 /2m 2 )). 

2.4.1 Relativistic Normalization 

We have defined the vacuum |0) which we normalize as (0| 0) = 1. The one-particle 




(2.51) 




(2.53) 



states \p) = at|0) then satisfy 



(p\q) = (27rr5^(p-q) 



(2.54) 
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But is this Lorentz invariant? It's not obvious because we only have 3-vectors. What 
could go wrong? Suppose we have a Lorentz transformation 

p" - (p'Y = KP U (2-55) 

such that the 3-vector transforms as p — > p' . In the quantum theory, it would be 
preferable if the two states are related by a unitary transformation, 

\p) - \p') = U(A) \p) (2.56) 

This would mean that the normalizations of \p) and \p') are the same whenever p and 
p' are related by a Lorentz transformation. But we haven't been at all careful with 
normalizations. In general, we could get 

\p)^X(p,p') \p') (2.57) 

for some unknown function \(p,p'). How do we figure this out? The trick is to look at 
an object which we know is Lorentz invariant. One such object is the identity operator 
on one-particle states (which is really the projection operator onto one-particle states). 
With the normalization (2.54) we know this is given by 

(0s Id (Pi ( 2 - 58 ) 

This operator is Lorentz invariant, but it consists of two terms: the measure J d 3 p and 
the projector \p)(p\- Are these individually Lorentz invariant? In fact the answer is no. 



Claim The Lorentz invariant measure is, 

<Pp 



r tfp_ 

J 2E? 



(2.59) 



Proof: J d A p is obviously Lorentz invariant. And the relativistic dispersion relation 
for a massive particle, 



p^ = m 2 pl = El = p 2 + m 2 (2.60) 

is also Lorentz invariant. Solving for p , there are two branches of solutions: po = ±Ef. 
But the choice of branch is another Lorentz invariant concept. So piecing everything 
together, the following combination must be Lorentz invariant, 

d 3 p 



[ d 4 p 5{p 2 - p 2 - m 2 ) = I 7T- 
J PO >o J z Po 



(2.61) 



which completes the proof. □ 



- 32 - 



From this result we can figure out everything else. For example, the Lorentz invariant 
5-function for 3-vectors is 

2E f 5 ( - 3) (p- q) (2.62) 

which follows because 

/^2%*< 3 >(jf-fl = l (2.63) 

So finally we learn that the relativistically normalized momentum states are given by 

^|p) = V^4|0) (2.64) 



v 

Notice that our notation is rather subtle: the relativistically normalized momentum 
state \p) differs from \p) just by the absence of a vector over p. These states now satisfy 

(p\q) = (2ixf2E l! 5^\p-q) (2.65) 

Finally, we can rewrite the identity on one-particle states as 

Some texts also define relativistically normalized creation operators by a\p) = ^/2Efal. 
We won't make use of this notation here. 

2.5 Complex Scalar Fields 

Consider a complex scalar field ip(x) with Lagrangian 

C = d^* - M 2 ^ (2.67) 

Notice that, in contrast to the Lagrangian (1.7) for a real scalar field, there is no factor 
of 1/2 in front of the Lagrangian for a complex scalar field. If we write ip in terms 
of real scalar fields by ip = (0i + i(j) 2 )/\/2, we get the factor of 1/2 coming from the 
l/\/2's. The equations of motion are 

ip + M 2 ip = 
d^i/j* + MV = (2.68) 

where the second equation is the complex conjugate of the first. We expand the complex 
field operator as a sum of plane waves as 
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Since the classical field ip is not real, the corresponding quantum field ip is not hermitian. 
This is the reason that we have different operators b and c' appearing in the positive 
and negative frequency parts. The classical field momentum is n = dC/dip = ip*. We 
also turn this into a quantum operator field which we write as, 

7T 

The commutation relations between fields and momenta are given by 

[ij(x),n(y)]= i 5^(x-y) and ^(x),n\y)] = (2.71) 

together with others related by complex conjugation, as well as the usual [ip(x),ip(y)] = 
[ip(x), ip'(y)] — 0, etc. One can easily check that these field commutation relations are 
equivalent to the commutation relations for the operators bp and cp, 

[b p ^ = (2n) 3 S^(p-q) 

[c^c^ = (2nf5^(p-q) (2.72) 

and 

[bp, bq] = [cp, c ? -] = [bp, c f ] = [bp, ct] = (2.73) 

In summary, quantizing a complex scalar field gives rise to two creation operators, &t 
and ct. These have the interpretation of creating two types of particle, both of mass M 
and both spin zero. They are interpreted as particles and anti-particles. In contrast, 
for a real scalar field there is only a single type of particle: for a real scalar field, the 
particle is its own antiparticle. 

Recall that the theory (2.67) has a classical conserved charge 

Q = i J d 3 x (ip*ip - ip*ip) = i J d 3 x (nip - ip*n*) (2.74) 
After normal ordering, this becomes the quantum operator 

Q = J (4<* - b f V) =N C -N b (2.75) 

so Q counts the number of anti-particles (created by c') minus the number of particles 
(created by tf). We have [H,Q] = 0, ensuring that Q is conserved quantity in the 
quantum theory. Of course, in our free field theory this isn't such a big deal because 
both iV c and iVj, are separately conserved. However, we'll soon see that in interacting 
theories Q survives as a conserved quantity, while N c and N b individually do not. 



d 3 p 
(27p 





t p -ip-x 



) 



(2.70) 
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2.6 The Heisenberg Picture 

Although we started with a Lorentz invariant Lagrangian, we slowly butchered it as we 
quantized, introducing a preferred time coordinate t. It's not at all obvious that the 
theory is still Lorentz invariant after quantization. For example, the operators <f>(x) 
depend on space, but not on time. Meanwhile, the one-particle states evolve in time 
by Schrodinger's equation, 

i*^ = H\ftt)) => W)) = e-^ t \p) (2.76) 

Things start to look better in the Heisenberg picture where time dependence is assigned 
to the operators O, 

O h = e iHt O s e- im (2.77) 

so that 

*^- = i[H,0 B ] (2.78) 

where the subscripts S and H tell us whether the operator is in the Schrodinger or 
Heisenberg picture. In field theory, we drop these subscripts and we will denote the 
picture by specifying whether the fields depend on space (f>(x) (the Schrodinger picture) 
or spacetime (f>(x,t) = <j)(x) (the Heisenberg picture). 

The operators in the two pictures agree at a fixed time, say, t — 0. The commutation 
relations (2.2) become equal time commutation relations in the Heisenberg picture, 

[0(f,t),0(y,t)] = [7r(f,t),7r(f,t)] = O 

y>(x,t),n(M]=iS^(x-if) (2.79) 

Now that the operator <f>(x) = (f>(x, t) depends on time, we can start to study how it 
evolves. For example, we have 

4> = i[H^\ = l -[f d'y^yf + V^yf + m 2 ^) 2 , <j>(x)} 

= i J d 3 y ir(y) (-i) 5 (3) (y - x) = ir(x) (2.80) 
Meanwhile, the equation of motion for n reads, 

* = i[H,*\ = \[J d 3 y ^y) 2 + V0(y) 2 + m 2 0(y) 2 , tt(x)] 
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+2im 2 (f)(y) 5 {3 \x-y) 

= V 2 0-m 2 (2.81) 

where we've included the subscript y on V y when there may be some confusion about 
which argument the derivative is acting on. To reach the last line, we've simply inte- 
grated by parts. Putting (2.80) and (2.81) together we find that the field operator <fi 
satisfies the Klein-Gordon equation 

d^(j) + m 2 (f) = (2.82) 

Things are beginning to look more relativistic. We can write the Fourier expansion of 
(f>(x) by using the definition (2.77) and noting, 

e im a § e~ im = a § and e iHt a\e~ im = e +iE ^ a* (2.83) 

which follows from the commutation relations [H, a-] = — Epa- and [H, at] = +Efa~. 
This then gives, 

*<*•'> = / 1& tM^^ 4'***) (2 ' 84) 

which looks very similar to the previous expansion (2.18) 
except that the exponent is now written in terms of 4- 
vectors, p ■ x = Eft — p • x. (Note also that a sign has 
flipped in the exponent due to our Minkowski metric con- 
traction). It's simple to check that (2.84) indeed satisfies 
the Klein-Gordon equation (2.82). 

2.6.1 Causality 

We're approaching something Lorentz invariant in the Figure 4: 

Heisenberg picture, where 4>(x) now satisfies the Klein- 
Gordon equation. But there's still a hint of non-Lorentz invariance because <fi and 7r 
satisfy equal time commutation relations, 

[<P(x,t),n(y,t)]=t5^(x-y) (2.85) 
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But what about arbitrary spacetime separations? In particular, for our theory to be 
causal, we must require that all spacelike separated operators commute, 



[O 1 (x),O 2 (y)] = V (x-y) 2 <0 (2.86) 

This ensures that a measurement at x cannot affect a measurement at y when x and y 
are not causally connected. Does our theory satisfy this crucial property? Let's define 

A(x-y) = [(j>(x),<l>(y)] (2.87) 

The objects on the right-hand side of this expression are operators. However, it's easy 
to check by direct substitution that the left-hand side is simply a c-number function 
with the integral expression 

A(x - y) = f (e- ip<x - y) ~ e*™) (2.88) 

What do we know about this function? 



• It's Lorentz invariant, thanks to the appearance of the Lorentz invariant measure 
j d 3 p/2E^ that we introduced in (2.59). 

• It doesn't vanish for timelike separation. For example, taking x — y = (t, 0, 0, 0) 
gives [<f>(x, 0), (j){x , *)] ~ e- imt - e +imt . 

• It vanishes for space-like separations. This follows by noting that A(x — y) = 
at equal times for all (x — y) 2 = — (x — y) 2 < 0, which we can see explicitly by 
writing 

W , t), m *)] = 5 / tSs t4t=2 ( eiP<S ~ Vl ~ ^ { ^) (2-89) 

and noticing that we can flip the sign of p in the last exponent as it is an inte- 
gration variable. But since A(x — y) Lorentz invariant, it can only depend on 
(x — y) 2 and must therefore vanish for all (x — y) 2 < 0. 

We therefore learn that our theory is indeed causal with commutators vanishing 
outside the lightcone. This property will continue to hold in interacting theories; indeed, 
it is usually given as one of the axioms of local quantum field theories. I should mention 
however that the fact that [(f>(x), <f){y)} is a c-number function, rather than an operator, 
is a property of free fields only. 
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2.7 Propagators 

We could ask a different question to probe the causal structure of the theory. Prepare 
a particle at spacetime point y. What is the amplitude to find it at point x? We can 
calculate this: 



(01^)10) = / ^ -^ (0 K4 ,|0) 



(2.90) 



The function D(x — y) is called the propagator. For spacelike separations, (x — y) 2 < 0, 
one can show that D(x — y) decays like 

D(x-y) ~ e - m \*-y\ (2.91) 

So it decays exponentially quickly outside the lightcone but, nonetheless, is non- vanishing! 
The quantum field appears to leak out of the lightcone. Yet we've just seen that space- 
like measurements commute and the theory is causal. How do we reconcile these two 
facts? We can rewrite the calculation (2.89) as 

[0(x), (f>{y)\ = D(x - y )-D(y-x) = if (x - yf < (2.92) 

There are words you can drape around this calculation. When (x — y) 2 < 0, there 
is no Lorentz invariant way to order events. If a particle can travel in a spacelike 
direction from x — > y, it can just as easily travel from y — > x. In any measurement, the 
amplitudes for these two events cancel. 

With a complex scalar field, it is more interesting. We can look at the equation 
[ip(x), ip^(y)] = outside the lightcone. The interpretation now is that the amplitude 
for the particle to propagate from x — > y cancels the amplitude for the antiparticle 
to travel from y — > x. In fact, this interpretation is also there for a real scalar field 
because the particle is its own antiparticle. 

2.7.1 The Feynman Propagator 

As we will see shortly, one of the most important quantities in interacting field theory 
is the Feynman propagator, 

A F (x-y) = (Q\T</,(x)<f,(y)\Q) = { ~ v) X ° > V [ (2.93) 

D(y — x) y v > x v 
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where T stands for time ordering, placing all operators evaluated at later times to the 
left so, 



T(f)(x)(p(y) 



x°>y° 
y°>x° 



(2.94) 



(f)(x)(f)(y) 
<j)(y)(j)(x) 

Claim: There is a useful way of writing the Feynman propagator in terms of a 4- 
momentum integral. 

d 4 p i 



A F (x - y) = 



■■(x-y) 



(2.95) 



(2tt) 4 p 2 

Notice that this is the first time in this course that we've integrated over 4-momentum. 
Until now, we integrated only over 3-momentum, with p° fixed by the mass-shell con- 
dition to be p° = Ep. In the expression (2.95) for A F , we have no such condition on 
p°. However, as it stands this integral is ill-defined because, for each value of p, the 
denominator p 2 — m 2 = (p ) 2 —p 2 —m 2 produces a pole when p° = ±Ep = ±^/p 2 + m 2 . 
We need a prescription for avoiding these singularities in the p integral. To get the 
Feynman propagator, we must choose the contour to be 

Im(po) 



Re(pO) 





Figure 5: The contour for the Feynman propagator. 

Proof: 

1 - I = I (2 96) 

p 2 -m 2 {p°) 2 -El (p° - E^){p° + E p ) V ' J 

so the residue of the pole at p° = ±Ep is ±l/2Ep. When x° > y°, we close the contour 
in the lower half plane, where p° — > —too, ensuring that the integrand vanishes since 
e -ip°(x°-y°) _^ q rj^g i^egral over p° then picks up the residue at p° = +Ep which 
is —2m/2Ep where the minus sign arises because we took a clockwise contour. Hence 
when x° > y° we have 



FK y> J (2tt) 4 2E p 



{3-y) 
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which is indeed the Feynman propagator for x° > y°. In contrast, when y° > x°, we 
close the contour in an anti-clockwise direction in the upper half plane to get, 

(2vr) 3 2Ej? 

0_A_ e - tr ^ = D(y _ x) (2 . 98) 

where to go to from the second line to the third, we have flipped the sign of p which 
is valid since we integrate over d 3 p and all other quantities depend only on p 2 . Once 
again we reproduce the Feynman propagator. □ 

Instead of specifying the contour, we may instead write im(p°) 
the Feynman propagator as 

A F (x -y)= -— — y— - 2.99 ^ 

J (27r) 4 p l — m z + ie I ~i£ 

with e > 0, and infinitesimal. This has the effect of shifting 

the poles slightly off the real axis, so the integral along the 

real p° axis is equivalent to the contour shown in Figure 5. Figure 6. 

This way of writing the propagator is, for obvious reasons, 

called the "ie prescription" . 

2.7.2 Green's Functions 

There is another avatar of the propagator: it is a Green's function for the Klein-Gordon 
operator. If we stay away from the singularities, we have 

(d? - V 2 + m 2 )A F ( X " V) = / (-P 2 + - 2 ) *- iHX - y) 



dp 

= —i I -, — rr e 



'4„. 

ip-(x-y) 



(2vr) 4 

= -i6W(x-y) (2.100) 

Note that we didn't make use of the contour anywhere in this derivation. For some 
purposes it is also useful to pick other contours which also give rise to Green's functions. 
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Im(pO) Im(p<>) 



Re(p °) ^ J ; ^_ > : 1^'<1' ,: ) 



Figure 7: The retarded contour Figure 8: The advanced contour 

For example, the retarded Green's function A R (x — y) is defined by the contour shown 
in Figure 7 which has the property 

f D(x -y)- D(y - x) x° > y° 
A R (x-y) = \ (2-101) 

[ y° > x° 

The retarded Green's function is useful in classical field theory if we know the initial 
value of some field configuration and want to figure out what it evolves into in the 
presence of a source, meaning that we want to know the solution to the inhomogeneous 
Klein-Gordon equation, 

d^(j) + m 2 (f) = J(x) (2.102) 

for some fixed background function J(x). Similarly, one can define the advanced Green's 
function Aa(x — y) which vanishes when y° < x°, which is useful if we know the end 
point of a field configuration and want to figure out where it came from. Given that 
next term's course is called "Advanced Quantum Field Theory", there is an obvious 
name for the current course. But it got shot down in the staff meeting. In the quantum 
theory, we will see that the Feynman Green's function is most relevant. 

2.8 Non-Relativistic Fields 

Let's return to our classical complex scalar field obeying the Klein-Gordon equation. 
We'll decompose the field as 

i/>(x,t) = e" im V(f,t) (2.103) 

Then the KG-equation reads 

dfy - V 2 ^ + m 2 ij = e~ irnt $ - 2iimp - V 2 ^ = (2.104) 

with the m 2 term cancelled by the time derivatives. The non-relativistic limit of a 
particle is \p\ m. Let's look at what this does to our field. After a Fourier transform, 
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this is equivalent to saying that <m|i|)|. In this limit, we drop the term with two 
time derivatives and the KG equation becomes, 

This looks very similar to the Schrodinger equation for a non-relativistic free particle of 
mass m. Except it doesn't have any probability interpretation — it's simply a classical 
field evolving through an equation that's first order in time derivatives. 

We wrote down a Lagrangian in section 1.1.2 which gives rise to field equations 
which are first order in time derivatives. In fact, we can derive this from the relativistic 
Lagrangian for a scalar field by again taking the limit d t ip <C rmjj. After losing the 
tilde, so ip — > ip, the non-relativistic Lagrangian becomes 

£ = - _ V^* W> (2.106) 

2m 

where we've divided by l/2m. This Lagrangian has a conserved current arising from 
the internal symmetry tp — > e ia ip. The current has time and space components 

f = (VV>, 2^(V>*W> - (2.107) 
To move to the Hamiltonian formalism we compute the momentum 

7l= ( ^l =l ^ (2.108) 

dip 

This means that the momentum conjugate to ip is iip*. The momentum does not depend 
on time derivatives at all! This looks a little disconcerting but it's fully consistent for a 
theory which is first order in time derivatives. In order to determine the full trajectory 
of the field, we need only specify ip and ip* at time t = 0: no time derivatives on the 
initial slice are required. 

Since the Lagrangian already contains a "pg" term (instead of the more familiar |pg 
term), the time derivatives drop out when we compute the Hamiltonian. We get, 

^=-Lv^V^ (2.109) 
2m 

To quantize we impose (in the Schrodinger picture) the canonical commutation relations 

[ij{x)^{y)}=5^{x-y) (2.110) 
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We may expand ip(x) as a Fourier transform 

m = I Wf a ^ (2 - m) 

where the commutation relations (2.110) require 

\a^a\} = {2^\p-q) (2.112) 

The vacuum satisfies a-|0) = 0, and the excitations are at^ . . . at |0). The one-particle 
states have energy 

H\P) = ^\P) (2.113) 

which is the non-relativistic dispersion relation. We conclude that quantizing the first 
order Lagrangian (2.106) gives rise to non-relativistic particles of mass m. Some com- 
ments: 

• We have a complex field but only a single type of particle. The anti-particle is 
not in the spectrum. The existence of anti-particles is a consequence of relativity. 

• A related fact is that the conserved charge Q = J d 3 x : ip'ip : is the particle 
number. This remains conserved even if we include interactions in the Lagrangian 
of the form 

AC = V{^) (2.114) 

So in non-relativistic theories, particle number is conserved. It is only with rela- 
tivity and the appearance of anti-particles, that particle number can change. 

• There is no non-relativistic limit of a real scalar field. In the relativistic theory, 
the particles are their own anti-particles, and there can be no way to construct a 
multi-particle theory that conserves particle number. 

2.8.1 Recovering Quantum Mechanics 

— * 

In quantum mechanics, we talk about the position and momentum operators X and 

— * 

P. In quantum field theory, position is relegated to a label. How do we get back to 
quantum mechanics? We already have the operator for the total momentum of the field 

p = jS? t ^ a * (2115) 
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which, on one-particle states, gives P \p) = p\p). It's also easy to construct the position 
operator. Let's work in the non-relativistic limit. Then the operator 

|iH e ~** (2,116) 

creates a particle with ^-function localization at x. We write \x) = ip^(x) |0). A natural 
position operator is then 

X = J d 3 x x ^(x) il){x) (2.117) 

so that X \x) = x\x). 

Let's now construct a state \ip) by taking superpositions of one-particle states \x), 

\ip) = I d 3 x <p(x) \x) (2.118) 



The function (p(x) is what we would usually call the Schrodinger wavefunction (in the 
position representation). Let's make sure that it indeed satisfies all the right properties. 
Firstly, it's clear that acting with the position operator X has the right action of <p(x), 

X 1 \<p) = J d 3 x x i (p(x) \x) (2.119) 

— * 

but what about the momentum operator PI We will now show that 

P>> = jd 3 x |f) (2.120) 

which tells us that P % acts as the familiar derivative on wavefunctions \(p). To see that 
this is the case, we write 

p * lv?) = / i0fH^v^)i>\t) |o) 

= / W^ e "" X>(f) |0) (2 ' 121) 

where we've used the relationship [a-,^(a?)] = e~ l ^' x which can be easily checked. 
Proceeding with our calculation, we have 

d?xcPp ifS ( d<p\ t 



( ( -i^- I a} |0) 



(2tt) 3 V dx i 

. dip 
dx i 



= (-'Si W (2-122) 
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which confirms (2.120). So we learn that when acting on one-particle states, the oper- 
ators X and P act as position and momentum operators in quantum mechanics, with 
[X\P^] \ip) = \<f). But what about dynamics? How does the wavefunction <p(x,t) 
change in time? The Hamiltonian (2.109) can be rewritten as 

H = f*z± W V* = / (gs £ «K ("23) 



so we find that 



.dtp 



2m 



vv 



(2.124) 



But this is the same equation obeyed by the original field (2.105)! Except this time, it 
really is the Schrodinger equation, complete with the usual probabilistic interpretation 
for the wavefunction tp. Note in particular that the conserved charge arising from the 
Noether current (2.107) is Q — J d 3 x \(f(x)\ 2 which is the total probability. 



Historically, the fact that the equation for the classical field (2.105) and the one- 
particle wavefunction (2.124) coincide caused some confusion. It was thought that 
perhaps we are quantizing the wavefunction itself and the resulting name "second quan- 
tization" is still sometimes used today to mean quantum field theory. It's important to 
stress that, despite the name, we're not quantizing anything twice! We simply quantize 
a classical field once. Nonetheless, in practice it's useful to know that if we treat the 
one-particle Schrodinger equation as the equation for a quantum field then it will give 
the correct generalization to multi-particle states. 



Interactions 

Often in quantum mechanics, we're interested in particles moving in some fixed back- 
ground potential V(x). This can be easily incorporated into field theory by working 
with a Lagrangian with explicit x dependence, 

£ = l ^"i V ^ V ^ _1/ ^^ (2.125) 

Note that this Lagrangian doesn't respect translational symmetry and we won't have 
the associated energy-momentum tensor. While such Lagrangians are useful in con- 
densed matter physics, we rarely (or never) come across them in high-energy physics, 
where all equations obey translational (and Lorentz) invariance. 
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One can also consider interactions between particles. Obviously these are only impor- 
tant for n particle states with n > 2. We therefore expect them to arise from additions 
to the Lagrangian of the form 

AC = xjj*(x) i/>(x) i/>(x) (2.126) 

which, in the quantum theory, is an operator which destroys two particles before creat- 
ing two new ones. Such terms in the Lagrangian will indeed lead to inter-particle forces, 
both in the non-relativistic and relativistic setting. In the next section we explore these 
types of interaction in detail for relativistic theories. 
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